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Abstract A periodic plasmonic meta-material was stud¬ 
ied using finite-difference time domain (FDTD) method 
to investigate the influence of neighboring particles on 
the near unity optical absorptivity. The meta-material 
was constructed as a silver nanoparticle (20-90nm) situ¬ 
ated above an alumina (AI2O3) dielectric environment. 
A full parametric sweep of the particle width and the 
dielectric thickness was conducted. Computational re¬ 
sults identified several resonances between the metal- 
dielectric and metal-air that have potential to broad¬ 
ening the response through stacked geometry. A sig¬ 
nificant coupled resonance between the metal-dielectric 
resonance and a cavity resonance between particles was 
capture as a function of dielectric thickness. This cou¬ 
pled resonance was not evident below dielectric thick¬ 
nesses of 40nm and above cavity widths of 20nm. Addi¬ 
tionally, a noticeable propagating surface plasmon po- 
lariton resonance was predicted when the particle width 
was half the unit cell length. 


1 Introduction 

Over the past decade there has been a tremendous 
interest in manipulating optical wavelengths by tailor¬ 
ing the typology of plasmonic meta-materials (MMs) [T]. 
The tailoring of electromagnetic MM’s have demon¬ 
strated augmented optical properties such as negative 
index of refraction, artificial permittivity, and perme¬ 
ability. More recently, research of meta-materials has 
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focused on two extremes of optical wavelength manipu¬ 
lation. These two extremes include, 1) the development 
of meta-materials with low-loss characteristics for op¬ 
tical switching applications m and 2) the develop¬ 
ment of meta-materials with nearly perfect absorption 
properties for energy conversion applications [illH] . The 
following research focuses on the latter by investigat¬ 
ing the influence of geometrical aspects on perfect ab¬ 
sorbers and understanding the neighboring interactions 
with future application for solar thermal energy conver¬ 
sion. 

One of the difficulties in understanding the response 
of a perfect optical absorbing meta-material is there are 
several geometry parameters that are critical to main¬ 
taining the desired response. Furthermore, as will be 
discussed in this study, as neighboring nanoparticles are 
brought close together, competing interactions act both 
destructively or constructively to either enhance or di¬ 
minish the response. By understanding the neighboring 
influences on the optical response there is a potential 
to broaden the perfect absorption characteristics. In the 
proposed application of this study, this equates to in¬ 
crease solar to thermal conversion efficiency. 

To limit the complexity of this study the plasmonic 
material geometry will be limited to a simple metallic 
nanoparticle situated on top of a dielectric thin film fol¬ 
lowed by a thick ground plane. This geometry has been 
demonstrated both experimentally and theoretically [SJ 
17115] to have a nearly perfect absorption response. Many 
of these authors also demonstrated that by controlling 
the particle diameter and dielectric thickness the re¬ 
sponse of these very simple devices can be tuned to 
particular optical frequencies EEHlin]. By taking the 
concept of the single particle perfect absorber and intro¬ 
ducing the neighboring particle it is hypothesized that 
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a third degree of freedom will allow the response to be 
broadened. 

Meta-materials have been both experimentally and 
computationally shown to respond to all wavelengths of 
electromagnetic radiation [12] but the visible region has 
garnered a lot of interest for its applicability to many 
optical applications. Because our objective is to design 
a solar absorber the frequency of interest will be limited 
to the visible regime, that is, from 428-750nm (400-700 
THz). Justification for focusing on this region initially 
can be reasoned based on a plot of the solar irradiance 
outside the earth’s atmosphere as show in Figure[T]and 
the corresponding left axis. The right axis of the same 
figure is the integrated spectral irradiance and can be 
more easily interpreted as the maximum heat flux if 
there is complete absorption of the incident electromag¬ 
netic radiation. The same figure also brings emphasis 
to the integrated irradiance of only the visible region, 
which is 560lT/m^. The visible region contributes to 
nearly 50% of the total heat flux when compared to the 
total heat flux of 1356.121 T/to^ for the whole spectrum. 

In comparison to conventional materials such as car¬ 
bon black and graphite, which can absorb heat in ex¬ 
cess of 85% [13] it is reasoned that a broadband meta¬ 
materials could potentially absorb heat in excess of 90% 
as demonstrated at single frequencies [HHllS]. Fur¬ 
thermore, a limitation of carbon based materials stems 
from the reflectivity at air-carbon interface and the 
resulting emissivity of approximately 0.8. Recent ap¬ 
proaches to use carbon nanotubes m have increased 
the emissivity upwards of 0.9 but there are practica¬ 
ble limitations to adhering nanotube to a substrates 
and providing good thermal transport across the in¬ 
terface. A unique advantage of the proposed metal- 
dielectric materials is the thermal barrier of the in¬ 
terface is avoided by exciting plasmons within the di¬ 
electric or alumina medium. Furthermore, alumina has 
an inherently large thermal conductivity (dObF/mAT at 
20C) due to the large phonon contribution to the ther¬ 
mal conductivity that will aid in the dissipation of lo¬ 
calized dielectric heating. 

In this study it is assumed that the incident elec¬ 
tromagnetic energy will be absorbed in the form of 
surface oscillating plasmons. Followed by the assump¬ 
tion that the plasmons will decay into lattice vibra¬ 
tion or phonons and diffuse through the substrate of 
the device. Often this decay of plasmons to phonons 
is referred to as dielectric heating but can also be re¬ 
lated to Joule heating. By assuming the material is 
dispersive and that the magnetic response is negligi¬ 
ble, it is possible to write a simple expression that 
relates the complex permittivity of a material to the 
volumetric energy generation. In order to utilize this 


Wavelength (nm) 



Frequency (THz) 


Fig. 1 Plot of the solar spectral irradiance as a function 
of frequency regenerated from data obtained from ASTM 
490 1161 . The right axis correspond to the integrated spectral 
irradiance, which has a cumulative total of 1356.12W/m^. 

expression, it is assumed that the material undergoes 
a forced harmonic oscillation with an incident electric 
field. Therefore, a material’s electronic response to a 
polarizing field can be characterized with the following, 
D = egE -t- P = egeE. This expression governs that the 
displacement (D) of the electrons is proportional to the 
applied electric field (E) and the polarization (P). The 
expression can be utilized to quantify the volumetric 
heating through Poynting’s theorem, which is derived 
in the literature [B]. Under the harmonic assumption, 
the time averaged expression for the volumetric heat 
generation within the material will take the following 
form. 


9 = iega;Im(e(a;))||E|p, (1) 

where eg is the permittivity of free space, u! is the 
angular frequency, ||E|| is the norm of the electric field 
vector, and Im(e(w)) is the imaginary part of the per¬ 
mittivity of the dispersive material. The units of this 
expression are in power (work done per unit time) per 
unit volume and the general form of the equation can 
be related to Joule heating (P(t) = I{t)V{t)). This pre¬ 
vious expression can be integrated over volume to de¬ 
termine the total power absorbed by the material at a 
particular frequency. 


q = j iegw/TO(e(a;))||Ef dU (2) 
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2 Plasmonic Material Design 

To limit the complexity of the antenna, a very sim¬ 
ple metal-dielectric-metal construction was implemented 
as illustrated in Figures [5J A unit cell representation 
was used to alleviate the computational expense of the 
simulation while still maintaining an accurate response 
of the whole material through the use of Bloch bound¬ 
ary conditions at the minimum and maximum x- and y- 
boundaries. The geometry of the nanoparticle anchored 
on top of a dielectric thin film was square with extent 
in X- and y-dir being equivalent to maintain the square 
geometry throughout the study. The thickness of the 
particle was held at 50nm in the z-dir. Below the dielec¬ 
tric material a thick ground plane (substrate) is situ¬ 
ated with the same composition of the top nanoparticle. 
The remaining space between the top of the material 
and the upper port was assumed vacuum region. 



Fig. 2 Dimensional illustration of the meta-material nano¬ 
antenna constructed from a nanoparticle separated from a 
ground plane with an alumina dielectric. Critical parame¬ 
ters to this research are the particle width {Wp) and dielec¬ 
tric thickness (To). The meta-material composite is perfectly 
symmetrical and infinite in the xy-plane. 


The top nanoparticle and the bottom ground plane 
for the following simulations were composed of silver 
(Ag). The dielectric thin film situated between the top 
particle and the bottom ground place was chosen to 
have a composition of alumina or aluminum oxide (AI 2 O 3 ). 
The substrate had dimensions Wq x Lq x Tq upon 
which the alumina with identical lateral dimensions Wq x 
WgxTjj was constructed. The silver nano-antenna with 
dimensions Wp x Lp xTp was located in the geometric 
center on top of the alumina layer. By controlling the 
extent of the substrate lateral dimension Wq x Lq the 
distance between neighboring particles was controlled. 

The boundaries that lie in the y-z and x-z plane 
were periodic such that the transverse magnetic field 
in the y-z plane was zero and the transverse electric 
field in the x-z plane was zero. A polarized TEM plane- 
wave was propagated in the —z direction from the open 
boundary in the x-y plane in the maximum z location 
of the computational domain. The bottom of the com¬ 
putational domain in the minimum z location had an 


open boundary condition. The scattering parameters 
were monitored at both of the boundaries in the z-dir. 

A Drude model was employed because the silver ma¬ 
terial is dispersive. Table [I] provides a list of Drude pa¬ 
rameters for silver. The alumina was assumed loss-less 
with a fixed permittivity of 6.20, which is equivalent 
to a refractive index of 2.49 (n = \/(e)). This value is 
based on experimental studies of alpha-phase alumina 
(a — AI 2 O 3 ) [17]. The permeability of all the materials 
was assumed to be unity provided this study is focused 
on the optical region of the electromagnetic spectrum. 

Because of the variation of the reported Drude pa¬ 
rameters in recent literature, the Drude parameters were 
compared to experimental data and a least squares fit 
was used to determine the Drude parameters for this 
study. In order to plot the Drude parameters the per¬ 
mittivity was written in the real and imaginary form. 
By assuming the material response to incident radia¬ 
tion is harmonic the Drude model that describes the 
complex permittivity takes the following form. 


e(w) = £00 --= ei + 162, (3) 

+ ILOJc 

where ujp is the plasma frequency, 7 c is the collision 
frequency. The complex permittivity can be separated 
into a real part (ei) and an imaginary part (£ 2 ). Fol¬ 
lowing the derivation of Ordal m. the real part and 
imaginary part take the following form. 


UJZ 


— ^00 




(27r7c)^ 


(4) 


a;(a;2 _|_ (27r7c)2) ■ 

Using the equations for ei and £2 the data provided 
in the experimental research of Christy |19j was fit with 
free parameter uip and 7 c. The result of this fit along 
with common parameters cited from other studies are 
depicted in Figure [3] and the corresponding Table [T| 
It is mentioned by Christy, there is a large error in 
£2 in the infrared region when the wave vector (fc) is 
large compared to n, stemming from the relation £2 = 
2nk. Christy goes on to mention that the error in £2 is 
less significant in the optical region where n and k are 
proportional. Therefore, for this study the experimental 
data is used with mild confidence but the reader should 
be aware of the apparent difference between £2 fit using 
the Drude model and the experimental data. 
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Source 

Plasma Freq. (ojp) [rad/s] 

Collision Freq. ( 7 c) [Hz] 

Baida et al. Loi 

1.37x10^^ 

S.lOxlQi-’ 

Ordal et al. 1181 

1.37x10^“ 

6.57x1011’ 

Fit to Christy et al. 1191 

1.37x101® 

4.97x101’! 


Table 1 Dispersion parameters for a Drude dispersion model to describe the dispersive nature of the silver material. All 
materials were assumed to have a permeability of unity. 


-100 


Wavelength (nm) 

1000 750 600 500 428 

' Real Epsifon (e-,) - Exp, Christy et al. 1972' +” 

Imag Epsilon (eg) - Exp, Christy et al. 1972 x 

Real Epsilon (e.,) - Drude, Baida et al. 2003 - 

Imag Epsilon (eg) - Drude, Baida et al. 2003 - 

Real Epsilon (e^) - Drude, Ordal et al. 1985 . 

Imag Epsilon (e?) - Drude, Ordal et al. 1985 -.. 

Rear Epsilon (e.,) - Least Squares Fit to Exp - 

Imag Epsilon (Sg) - Least Squares Fit to Exp - 





+ -I- ^ 


Binf - 1.0 

cOp fit = 1.37E-I-16 rad/s 
7 cfit = 4.97E+12Hz 


3.5 
3 

2.5 


400 500 600 

Frequency (THz) 


Fig. 3 Plot of the real and imaginary permittivity comparing 
the Drude model to experimental data m- The inset param¬ 
eter (ojpfit and 7 cfit) were determined through a least squares 
fit of Christy’s experimental data. 


3 Computational Model 


A finite-difference time-domain (FDTD) computa¬ 
tional method employed in the CST Microwave Studio 
Package was used to predict the scattering parameters 
as a function of neighboring interactions and geome¬ 
try of the antenna. To aid in relating the scattering 
parameters to physically measurable quantities the sil¬ 
ver substrate of the device was constructed sufficiently 
thick to eliminate the transmission of electromagnetic 
energy. 

The absorptivity and the reflectivity from the view 
point of the top most boundary (port) in the positive 
z-dir above the device were predicted through the scat¬ 
tering parameters at that port. Because the silver sub¬ 
strate had a sufficient thickness the scattering param¬ 
eter from the top port to the bottom port (S 12 ) was 
assumed and confirmed by the simulation to be zero. 
This resulted in simplification of the scattering param¬ 
eters’ relation, S'!! + S 12 = 1. Defining an energy bal¬ 
ance where, A(w) -|- R(a;) -I- T(a;) = 1, where A(a;) 
is the absorptivity, R(a;) is the reflectivity and T(a;) 
is the transmissivity, and assuming the transmissivity 
is zero the energy balance simplified to the following, 
A(a;) + R(a;) = 1. The reflectivity then can be defined 


as scattering parameter R{uj) = |S'ii|^ and hence the 
absorptivity can be defined as A(a;) =l-|S'ii|^. 

A tetrahedral mesh was used to discretize the com¬ 
putational domain. The FDTD solution was specified 
to have ten grid points per wavelength. Over the fre¬ 
quency range up to 100 points were simulated and based 
on those points a 1001 point spline was fit to provide 
a smooth response over the range of interest. Also, the 
adaptive mesh refinement feature was used in CST to 
locally refine mesh cells and check grid mesh indepen¬ 
dence from the solution. 


4 Discussion 

This research was divided into two areas of interest, 
1 ) study of the optimal dielectric thickness and parti¬ 
cle diameter for a fixed domain of 91nm and 2) study 
of the influence of neighboring particles with varying 
domain length. The following section will begin with 
the optimization study of particle width and dielectric 
thickness. 


4.1 Absorptivity as a Function of Particle Width and 
Dielectric Thickness 

The absorptivity as a function of particle geometry 
and dielectric thickness was investigated by conducting 
a full parametric design sweep of the space. The ap¬ 
proach taken was to use the FDTD model to evaluate 
a series of 10,000 design points with the space to de¬ 
termine the trends of the design space to guide in the 
design of meta-material with peak absorptivity. This 
parametric study was carried out by fixing the domain 
width to 91nm. The first degree of freedom in the para¬ 
metric sweep was the particle width (Wp), which was 
varied from 20nm to 90nm for 100 samples. The sec¬ 
ond degree of freedom was the dielectric thickness (T^), 
which was varied from 20nm to lOOnm for 100 samples. 
The absorptivity was calculated for each configuration 
point and the results are shown in Figured which illus¬ 
trate the response of the peak absorption as a function 
of the two degrees of freedom. 

Figured] provide a complete view of the particle and 
dielectric response for the simple antenna. Captured in 
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Absorptivity 

o;o 



Fig. 4 Plot of the absorptivity as a function of particle width 
(Wp) and the dielectric thickness (T^) at a fixed domain 
length of 91nm (Ws)- Only points with an absorptivity above 
0.5 are shown in the figure. Dashed line 1 is the response due 
to the metal-air. Dashed line 2 is a coupled response of the 
metal-dielectric and the neighboring cavity resonance. At a 
dielectric thickness of 90nm the dashed line approaches the 
LSSP peak of the metal-dielectric. Cross-sections from this 
3D plot can be found in Figure [7] at a constant dielectric 
thickness and Figure [5] at a constant particle width. 

this figure is the response of not only the localized sur¬ 
face plasmon polaritons (LSPPs) due to both the metal- 
dielectric response and the metal-air response but also 
propagating SPPs (PSPPs). There are two trends that 
are most obvious and they lie on either plane of particle 
width V. frequency and dielectric thickness v. frequency. 
The first trend is outlined by the dashed line and the 
number 1 in the particle width v. frequency plane. This 
response is independent of the particle width and in¬ 
vestigating the spatial electric field plot determined the 
FDTD simulation the field is localized at the top cor¬ 
ners of the particles suggesting this a response due to 
the metal-air interface, where we assumed air has an 
electric permittivity of unity. The second trend on the 
adjoining plane of dielectric thickness v. frequency is 
a more interesting trend as it pertains to the metal- 
dielectric response. 

To aid in the interpretation of Figure 0] two differ¬ 
ent 2D slices have been taken. The first slice is in a 
dielectric thickness v. frequency plane at a fixed dielec¬ 
tric thickness of 23nm, see Figure [S] Within this slice 
there are three distinct region as labeled in the figure. 
The first region looking from left to right in the figure 
is the LSPPs that are at the four corners of the particle 
at the interface of the metal and dielectric (LSSP-MD). 
Following the 40nm particle width line, emphasized as 


a white dashed line in the figure the left most branch 
(416THz) is associated with a counter-clockwise circu¬ 
lar current ninety degrees out of phase with the inci¬ 
dent signal. The neighboring LSSP branch (497THz) 
is a clockwise circular current ninety degrees out of 
phase with the incident signal. In both of these cases 
Tithe heating is subsequently localized at the corners of 
CD the nanoparticle within the dielectric. Following further 
right to 582THz there is a propagating SPP resonance 
^ that propagates ninety degrees out of phase along the 
^ y-dir. The heating at this resonance would be along the 
..^x-min and x-max edge of the nanoparticle at the metal 
^dielectric interface. This is an advantageous resonance 
N from a thermal point of view as there would be an in¬ 
duced thermal gradient within the dielectric. The final 
area of interest along the line at (3) is the resonance 
due to the metal-air. All three of these resonances are 
visualized by electric field vectors in three dimensions 
in Figure El 

Focusing on the first point of interest in Figure El 
the dependence of the LSSP on the particle width pro¬ 
vides evidence for a stacked design that would result 
in broadband absorptions. If you were to form a pyra¬ 
midal stack m of alternating metal-dielectric layers 
where the dielectric was maintained a 23nm and the 
metal width was varied from 20nm to 60nm, the fre¬ 
quency range between 400THz and 550 THz could be 
absorbed. More interesting, there is a noted blue shift 
in the PSSP-MD resonance (Number 2 in Figure El as 
a function of dielectric thickness, which suggests that 
several 40nm particles atop increasing thickness of di¬ 
electric (20-90nm) would encompass a frequency range 
between 575THz and 4I6THz. Furthermore, inducing a 
thermal gradient across the plasmonic device. 

The second slice of interest from the 3D plot in Fig¬ 
ure 01 is in a dielectric thickness v. frequency slice at a 
particle width of 84nm, see Figure [71 Here we gather 
another range of resonances similar to what was shown 
in the previous slices. However, a strong dependence 
on dielectric thickness is noted and emphasized by the 
dashed line. Again, if we concentrate our discussion 
to a dielectric thickness of 40nm and start from the 
left and increase in frequency the first peak, which is 
not very significant at this geometry, is the LSSP res¬ 
onance at 4I6THz. This LSSP resonance is not signif¬ 
icant in intensity and lacks a dependence on the di¬ 
electric thickness. If a vector plot of electric field were 
investigated there would be a localized field at the cor¬ 
ners of the nanoparticle. Moving right, at (2) the reso¬ 
nance (58ITHz) is an interaction between the neighbor¬ 
ing particles within the cavity, which is often referred 
to as a cavity resonance. (3) is a PSSPs resonance that 
is influenced by the dielectric thickness. (4) is a LSSP- 
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Wavelength (nm) 

1000 856 750 666 600 545 500 461 428 400 375 



200 300 400 500 600 700 800 

Frequency (TFIz) 


Wavelength (nm) 

1000 856 750 666 600 545 500 461 428 400 375 



200 300 400 500 600 700 800 

Frequency (TFIz) 


Fig. 5 Plot of the (a) reflectivity and (b) absorptivity as a 
function of particle width and frequency. The dielectric thick¬ 
ness is T;;=23nm. Number (1) indicates localized surface plas- 
mon polaritons between the metal and dielectric (LSSP-MD). 
Number (2) is a resonance of propagating surface plasmon 
polaritons at the metal-dielectric interface (PSSP-MD). This 
PSSP-MD is peak when the particle width is half the domain 
length. Number (3) is a resonance of localized SSP at the 
metal-air interface (LSSP-MA) 


MA at the metal-air interface and (5) is a second or¬ 
der cavity resonance with the neighboring particles. If 
the dielectric thickness is greater than 40nm there is a 
strong interference that leads to the fundamental LSSP 
resonance at 416THz. In reference to Figure 0] this in¬ 
terference is a result of the response of the nanoparticle 
and the ground plane and the nanoparticle and neigh¬ 
boring nanoparticles. As the particle width decreases 
the interference is not apparent and the fundamental 
LSSP resonance is only present. 


4.2 Influence of Neighboring Particles on the 
Electromagnetic Response 

One of the complications with designing a periodic 
meta-material is the influence of neighboring structures. 
As noted in the previous section, the response of the 
material can be dependent on distance between neigh¬ 
boring particles. For this sub-study the influence of 
neighboring particles is discussed in terms of several 
known theories. In doing so, a series of simulations were 



Fig. 6 Plot of the electric field at three points of frequency 
that correspond to the three areas of emphasis in Figure [h] 
(a) is the LSSP-MD at 416THz, (b) is PSSP-MD at 582THz 
and, (c) is LSSP-MA at 714THz. 


conducted by increasing the distance between the edge 
of the particle and the periodic boundary condition for 
a single unit cell. The domain was fixed at 500nm and 
carried out in an identical manner as discussed in Sec¬ 
tion 01 

A series of four designs illustrated in Figure [8] were 
investigated. The first design involves varying both the 
particle spacing and dielectric thickness and defines the 
fundamental plasmon frequency. The geometries where 
the particle spacing and dielectric thickness have nearly 
perfect absorption is explained by anti-parallel current m- 
Figure 0] shows the absorption parameters for a select 
number of cases. As the period size decreases the funda¬ 
mental plasmonic frequency red shifts. This will serve 
as a baseline to identify the fundamental mode from 
the coming designs. 

The second design involves forming a cavity between 
the neighboring particles. This corresponds to Design B 
in Figure [HI The cavity response of this design is anal¬ 
ogous to that of a nanohole array where the incident 
electromagnetic wave resonates between the particles. 
Because of the proximity of the closely spaced silver 
nano-particles, extreme coupling occurs. This coupling 
leads to strong field enhancement in the near field, thus 
enhancing the light-matter interaction and resulting in 
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Fig. 7 Plot of the (a) reflectivity and (b) absorptivity as a 
function of dielectric thickness and frequency. The particle 
width is Wp=84nm. Number (1) is a weak resonance of the 
fundamental LSSP between the metal and dielectric. Number 
(2) is a cavity resonance, (3) is a PSSP resonance, (4) is a 
LSSP between the metal and air, and (5) is a second order 
cavity resonance. The angled dashed line above a dielectric 
thickness of 40nm is a result of coupling between the cavity 
resonance and metal-dielectric resonance. 



Fig. 8 Illustration of the four dominate interactions that are 
responsible for the response. Design A is an anti-parallel cur¬ 
rent design. Design B is a cavity interaction between particles. 
Design C is an effective medium design. Design D is a grating 
coupling design. 


perfect absorption with multiple resonances. Figure [TU] 
shows select runs where the distance between the par¬ 
ticles is decreased. As the distance decreases, the re¬ 
sponse red shifts. However, the fundamental mode re¬ 
mains stationary and the secondary mode is the peak 
that is shifted. This can be conceptually thought of 



Frequency (THz) 


Fig. 9 Results from Design A of an anti-parallel perfect ab¬ 
sorber. As the distance between particles increases the re¬ 
sponse red shifts. 


as longer wavelengths coupling within larger hole di¬ 
ameters. Additionally, the dielectric thickness does not 
influence the response as illustrated by the selected 
points. 



Frequency (THz) 

Fig. 10 Results from Design B of a cavity response between 
particles resulting in perfect absorption. As the cavity in¬ 
creases the response blue shifts. 


The third design is illustrated in Figure[5](Design 3). 
It is based on an effective medium theory. This theory 
presumes that the particles situated on top of the di¬ 
electric act as a singular material with material proper¬ 
ties that depend on the volume fraction of silver and air. 
As the distance between the particles increases (particle 
size decreases) the volume fraction of air increases. Per¬ 
fect absorption can be obtained by impedance matching 
the effective medium with the dielectric. A resonance of 
768.5 THz occurs when the period is 85nm (a lOnm dis¬ 
tance between particles) and the dielectric thickness is 
greater than 19nm. 
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To impedance match the effective medium with the 
dielectric, the ground plane was removed and a set par¬ 
ticle geometry was simulated. The peak frequency with 
a reflection of zero marked the impedance matching at 
768.5THz characterized by Z = \/^, where equals 
unity. Once this peak impedance-match frequency was 
known, the ground plane was put back in the simula¬ 
tion. The width of the domain was restricted to 85nm 
and the dielectric thickness was varied until the equiv¬ 
alent peak was found. The arrangement that demon¬ 
strated this similar response with the ground plane was 
a particle width of 85nm. This configuration resulted 
in the correct volume fraction of air and silver for the 
corresponding permittivity of air. A series of runs il¬ 
lustrate this in Figure [TTJ The response of this effective 
medium requires that the dielectric thickness be greater 
than 19nm to facilitate a long enough traverse time to 
absorb the energy. 



Frequency (THz) 


Fig. 11 Results of Design C of an effective medium material 
where the volume fraction of the silver and air impedance 
match the surrounding air. 

The fourth design is a grating coupling design that 
is illustrated by Design 4 in Figure [5] The grating cou¬ 
pling is evident by removing the ground plane and still 
maintaining the optical response. The grating coupling 
is realized when the period of between the particles is 
twice the particle width. This is illustrated in Figured^ 
where the particle width is set to 75nm and the particle 
period is ISOnm. The result of the ground plane with 
the grating results in a perfect absorber. 


5 Conclusion 

By understanding the interaction between neighbor¬ 
ing particles there is potential to design more efficient 
plasmonic devices but also the potential to broaden 



Frequency (THz) 


Fig. 12 Results of Design D of a grating coupler where the 
peak resonance is apparent when the spacing between the 
particles is twice the particle diameter. 

their response. In conducting a full parametric sweep of 
particle diameter and dielectric thickness at a fixed do¬ 
main length of 91nm there were several resonances that 
could be attributed to the metal-dielectric and metal¬ 
air response. At particle diameters within 20nm of the 
unit cell length, a coupled interaction between the cav¬ 
ity resonance and metal-dielectric resonance was noted. 
Similarly, for sufficiently thin dielectrics there was a de¬ 
pendence of the localized surface plasmon polaritons on 
the particle width at the corners of the nanoparticle 
and a strong resonance between the metal and air. The 
interactions of neighboring particles were reasoned to 
be explained by four current theories, which include 1) 
anti-parallel currents, 2) cavity resonances, 3) effective 
medium, and 4) grating coupling. Given the identifi¬ 
cation of these four interactions the influence on their 
geometry identified how the peak resonances shift. 
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